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The leading particle effect in charm hadroproduction is an enhancement of the cross section for 
a charmed hadron D in the forward direction of the beam when the beam hadron has a valence 
parton in common with the D. The large D + /D~ asymmetry observed by the E791 experiment is an 
example of this phenomenon. We show that the heavy-quark recombination mechanism provides an 
economical explanation for this effect. In particular, the D + / D~ asymmetry can be fit reasonably 
well using a single parameter whose value is consistent with a recent determination from charm 
photoproduction. 



Fixed-target hadroproduction experiments have ob- 
served large asymmetries in the production of charmed 
mesons and baryons These asymmetries are com- 

monly known as the "leading particle effect" , since 
charmed hadrons having a valence parton in common 
with the beam hadron are produced in greater numbers 
than other charmed hadrons in the forward region of the 
beam. For example, the E791 experiment Q, in which a 
500 GeV TT~(ud) beam is incident on a nuclear target, ob- 
serves a substantial excess of D~(cd) over D + (cd) when 
the charmed mesons are produced with large momentum 
along the direction of the tt~ beam. The asymmetry, 



da[D~] - da[D+] 
da[D~] + da[D+] 



(1) 



is as large as ~ 0.7 for the most forward D mesons 
measured. Asymmetries in the production of charmed 
baryons have also been observed |2| []|] . 

In contrast with the large experimental asymmetries, 
perturbative QCD predicts that the asymmetries should 
be very small. The QCD factorization theorem for heavy 
particle production || states that the cross section for 
producing a D meson in the collision of two hadrons can 
be written as 



da[hti -» D + X] 



(2) 
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where is the distribution function for the parton i in 
the hadron h, da(ij — * cc + X) is the parton cross section 
and D c ^d is the fragmentation function for a c quark 
hadronizing into a D meson. The corrections to Eq. (Q) 
are suppressed by powers of Aqcd/wc or Aqcd/p_l- A 
study of 0(Aqcd/"t-c) corrections to charm production 
can be found in Ref. || . The leading order (LO) diagrams 
for gg — > cc and qq — > cc produce c and c quarks symmet- 
rically. The c and c fragment independently into D and D 



mesons, and charge conjugation invariance requires that 
D c ^d = D-_^. Therefore, LO perturbative QCD pre- 
dicts no asymmetry between D and D mesons. Asymme- 
tries are generated at next-to-leading order (NLO) [f}-{l0) , 
but they are too small by an order of magnitude or more 
to account for the experimentally observed asymmetries. 

Thus, charm asymmetries are interesting because they 
probe the power corrections to Eq. (||). Most attempts 
to explain the leading particle effect have employed phe- 
nomenological models of hadronization. One explana- 
tion is the "beam drag effect" jll]] in which charm 
quarks hadronize into charmed mesons by the decay of a 
color string connected to partons in the beam remnant. 
Calculations of the beam drag effect employ the Lund 
string fragmentation model p"2| , as implemented in the 
PYTHIA Monte Carlo Q . Predictions for the asymme- 
tries are sensitive to the model for the beam remnant, 
as well as other parameters in PYTHIA. These models 
can be tuned to fit data if one uses a large charm quark 
mass and gives the partons intrinsic transverse momen- 
tum of ~ 1 GeV |lj . Another class of models |l4| are 
based on the possibility of intrinsic charm in the nucleon 
|l5| . These models are sensitive to the poorly determined 
intrinsic charm structure function as well as the model of 
the beam remnant and predict smaller asymmetries than 
are experimentally observed. 

In this letter we show that the leading particle ef- 
fect can be explained quantitatively and economically by 
the heavy-quark recombination mechanism introduced in 
Ref. This mechanism has already been applied to 

charm asymmetries in fixed-target photoproduction ex- 
periments in Ref. fL7| , where it was shown that the charm 
asymmetries observed in the experiments E687 and E691 
at Fermilab can be fit well with just one free parameter. 
Here we apply the heavy-quark recombination mecha- 
nism to the much larger asymmetries observed in the 
E791 hadroproduction experiment. 
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FIG. 1. Example of a diagram for cq recombination into a 
D meson. Single lines are light quarks, double lines are heavy 
quarks and the shaded blob is the D meson. 



In the heavy-quark recombination mechanism, a heavy 
quark or antiquark recombines with a light parton that 
participates in the hard-scattering process rather than a 
spectator parton from the beam or target. A nonpertur- 
bative parameter that characterizes the probability for 
the light parton to bind to the heavy quark to form a 
meson sets the overall normalization of the cross section. 
Since the asymmetry is generated in the hard process, 
the dependence of the asymmetry on kinematic variables 
is calculable within perturbative QCD. Therefore, heavy- 
quark recombination provides a simpler and more predic- 
tive explanation for the asymmetries than conventional 
hadronization models. 

An example of a Feynman diagram contributing to this 
process is depicted in Fig. [I]. A light antiquark q partic- 
ipates in a hard-scattering process which produces a c 
and c. The outgoing q emerges from the hard-scattering 
process with momentum that is soft in the rest frame of 
the outgoing c quark, and the c and q then bind to form a 
D meson. There is an analogous process in which a light 
quark combines with a c to form a D meson. We em- 
phasize that heavy-quark recombination is not taken into 
account by higher order perturbative QCD corrections to 
the fragmentation contribution. The NLO correction in- 
cludes a parton subprocess qg — > qcc that is similar to the 
parton subprocess of the diagram in Fig. [l]. However, for 
most of the phase space the outgoing q hadronizes into 
a jet that is distinct from the jets containing the c and 
c. When the q has small momentum in the c rest frame, 
nonperturbative effects can bind the c and q and enhance 
the cross section. This enhancement, which is not present 
in the NLO perturbative calculation, is accounted for by 
heavy-quark recombination. 

The heavy-quark recombination contribution to the D 
meson cross section is: 



da[D}=da[qg^(cq) n + c}p[(cqy 



D] 



(3) 



This cross section must be convolved with parton distri- 
butions for the q and g in the colliding hadrons. The sym- 
bol (cq) n indicates that the q has momentum O(Aqcd) 
in the rest frame of the c and that the (cq) n has the color 



and angular momentum quantum numbers specified by 
n. da(qg — ► {cq) n + c) is calculable in perturbative QCD, 
while p[(cq) n — > D] is a nonperturbative factor propor- 
tional to the probability for the {cq) n to evolve into a 
state including the D meson. 

A detailed description of the calculation of da\qg — » 
(cq) n + c] can be found in Ref. fujfl . In this letter, we 
summarize the most important features. In most regions 
of phase space, the heavy quark recombination contri- 
bution is power-suppressed relative to the fragmentation 
contribution in Eq. (|^) in accord with the factorization 
theorems of QCD. Let 9 be the angle between the incom- 
ing q and the outgoing D meson in the parton ccnter- 
of-momentum frame. At 9 — 7r/2, the ratio of the par- 
ton cross sections for heavy-quark recombination and LO 
gluon-gluon fusion is 



da[qg ^ (cq^S^) +c] 



da[gg -> cc] 
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Using the expected scaling p ~ Aqcd/w c , we find that 
heavy-quark recombination is an 0(AqcdWc/ 'p±) power 
correction for p± m c . However, the heavy-quark re- 
combination cross section is strongly peaked when the D 
is produced in the forward region of the incoming q. At 
9 = 0, the ratio of the parton cross sections is 



dafig^icq)^^^) 



da[gg 



8=0 
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81 
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Thus, there is no kinematic suppression factor when the 
(cq) n is produced in the forward region of the incoming q. 
Heavy-quark recombination provides a natural explana- 
tion for the leading particle effect because in the forward 
region, the cross section for a D{D) meson carrying a 
q{q) that matches a valence parton in the beam is larger 
than that of other charmed mesons since valence quark 
structure functions are larger than sea quark structure 
functions. 

If the (cq) n hadronizes into a D meson and nothing 
else, the (cq) n must be in a color-singlet state with the 
same flavor and angular momentum quantum numbers 
as the D. This contribution to p[(cq) n — > D\ is pro- 
portional to the square of a moment of the D meson 
light-cone wavefunction and it scales as Aqcd/wc- How- 
ever one should also allow for nonperturbative transi- 
tions in which the (cq) n hadronizes into states includ- 
ing additional soft hadrons in the rest frame of the D. 
The inclusive parameter p[{cq) n — > D], which includes 
these additional contributions, is still expected to scale 
like Aqcd/^c- Furthermore, the (cq) n can have differ- 
ent color and angular momentum quantum numbers than 
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the D meson in the final state. Since the momentum of 
the light quark is O(Aqcd), amplitudes for production 
of (cq) n with L > are suppressed relative to S-waves 
by powers of Aqcd/ to c or ^qcb/p±- The heavy-quark 
recombination contribution to D meson production thus 
contains four free parameters: 

pT = picdCS^) -» D+],pf = pM( 3 5«) - D+], (4) 
P r = P [cdCS ( s) ) - D+],pf = p[cd( 3 s{ 8) ) -> 

Here the superscript sm stands for spin-matched and sf 
stands for spin-flipped. (Nonperturbative transitions in 
which the light quark flavor of the (cq) n differs from that 
of the D are suppressed in the large N c limit of QCD and 
will be ignored in this analysis.) Heavy-quark symmetry 
relates the four parameters appearing in the production 
of D* + to the four parameters in Eq. (^). Parameters for 
D~ and D° production can be related using charge con- 
jugation and isospin symmetry. The color-octet contri- 
butions are analogous to the color-octet production ma- 
trix elements that play a prominent role in the theory 
of quarkonium production |l8| . If the d appearing in 
Eq. (Q) were a heavy antiquark, the parameters pf , p| m 
and pg would be suppressed by powers of v, where v is 
the typical velocity of the antiquark in the bound state. 
However, the d is light, so v ~ 1 and there is no apparent 
suppression of these parameters. 

In photoproduction, the 0(a a 2 a ) color-octet and color- 
singlet heavy-quark recombination cross sections have 
the same functional form, so the cross section depends 
only on linear combinations of color-singlet and color- 
octet parameters of the form pi + pg/8. The best fit 
to fixed-target photoproduction asymmetries in Ref. ]T^ ] 
yielded p\ m + /o| m /8 = 0.15. Including the p sf parame- 
ters does not improve the fit unless pf + p$ /8 is negative, 
which is unphysical. In the case of fixed target hadropro- 
duction, the color-octet and color-singlet cross sections 
have different functional forms, so the asymmetries de- 
pend on all four parameters. However, we find that the 
asymmetries measured by the E791 experiment can be fit 
rather well by the single parameter p\ m . In this analysis, 
we only include p\ m and set the other parameters to zero. 

There are two contributions to the heavy-quark recom- 
bination cross section for a D meson: 

a) ^ dcj[qg — > (eg)" + c] p[(cq) n — > D], (5) 

n 

b) J2 dd ^9 -> (cq) n + c] p n [(cq) n -^D]® D C ^ D . 

In process a), a light antiquark recombines with a c to 
form a D meson. In process ft), ac participates in the 
recombination process, and a D meson is produced via 
fragmentation from the c quark which does not recom- 
bine. For D production, process b) partially dilutes 
the asymmetry generated by a). We take m c — 1.5 GeV 
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FIG. 2. The asymmetry variable a[Z) + ] vs. xf and p\. 
The data points are the measurements of E791 and the solid 
line is our prediction with pi™ = 0.06. 

and set the renormalization and factorization scales to 
be \/p]_ + m%. The parton distributions are GRV 98 LO 
H for the nucleon and GRV-P LO @ for the pion. The 
E791 experiment uses a target consisting mostly of car- 
bon, so the nucleon structure function is an isospin sym- 
metric linear combination of proton and neutron. We use 
the one-loop expression for a s with 4 active flavors and 
Aqcd = 200 MeV. We include only the LO fragmenta- 
tion diagrams. If the NLO corrections are approximated 
by including a K factor, then pf m should be multiplied by 
the same K factor to obtain the same asymmetry. We use 
the Petersen parametrization for D c ^r> with e = 0.06 and 
the fragmentation probabilities are determined from data 
in Ref. Contributions to the D cross section from 

feeddown from D* decay are included. Feeddown from 
other excited D meson states is expected to be small. 

The predictions for a[D + ] as a function of xf and p\ 
are compared to the E791 data in Fig. 0. A least squares 
fit to all the data yields p| m = 0.06. The p± dependence 
is described very well by the heavy-quark recombination 
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mechanism. The xp distribution is described well for 
xp > 0.2, while our calculation underpredicts the asym- 
metry in the region 0.0 < xp < 0.2. It is possible to 
obtain a better fit to the asymmetry in this region if 
one allows nonvanishing values of pf , p| m and pf . While 
we expect heavy-quark recombination to dominate in the 
forward region, other 0(AQCD/m c ) corrections to frag- 
mentation may be important near xp — and could also 
account for the discrepancy. A systematic analysis of all 
0(Aqcd/"t-c) corrections to D cross sections is needed to 
address this issue, but this is beyond the scope of this pa- 
per. Although a was not measured at large negative xp, 
it must be positive in this region in accord with our pre- 
diction because of the leading particle effect associated 
with the target nucleons. 

The value of p\ m extracted from the fits is consistent 
with the scaling p\ m ~ 0(Aqcd/?71c)- A lower bound on 
pf m can be obtained by assuming that Pi m [(cq) n — > D] is 
saturated by a final state consisting of only the D. Then 
p\ m is proportional to the square of a moment of the 
D light-cone wavefunction which can be bounded using 
heavy-quark effective theory arguments p3| . This value 
is smaller than the lower bound pf m = 0.15 from pho- 
toproduction fits ]l7j] . However, the uncertainty in this 
extraction due to higher order perturbative corrections 
and 0(Aqcd/™c) nonperturbative corrections could eas- 
ily be a factor of 2. Furthermore, the fits in Ref. jl7|] used 
E687 data that was not corrected for efficiencies because 
of correlations between the measured asymmetries and 
the production model used in the Monte Carlo for simu- 
lating the trigger p3]. Given all these uncertainties, we 
regard the values of p| m extracted from hadroproduction 
and photoproduction to be consistent. 

It would be interesting to extend the analysis of this 
paper to asymmetries measured in experiments with dif- 
ferent beams and for different charmed hadrons, partic- 
ularly charmed baryons. Baryon asymmetries could be 
generated by a process similar to that of Fig. [I] with the 
light antiquark replaced by a light quark. There is an at- 
tractive force between the c and q when their color state 
is the 3 of SU(3) that can enhance the cross section in 
this channel. It is natural to expect the (cq) n diquark 
to hadronize to a state with a charmed baryon. The 
SELEX collaboration has measured Af asymmetries in 
experiments with 7r~, E~ and p beams [||. The A+ and 
A~ both have a valence quark that matches one of the 
valence quarks of the tt~ , while only the A+ has a valence 
quark in common with the p and £~. Therefore, heavy- 
quark recombination predicts that A[r asymmetries are 
much greater for p and S _ beams than for a tt~ beam. 
This prediction is in agreement with SELEX measure- 
ments, which find a:[A~] rs 1 for p and S _ beams and 
a[A~] 0.2 for ir~ beams. 
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